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CO ' A curious feature of complex scattering potentials v(x) is the appearance of spectral sin- 

gularities. We offer a quantitative description of spectral singularities that identifies them 
0^ \ with an obstruction to the existence of a complete biorthonormal system consisting of the 

eigenf unctions of the Hamiltonian operator, i.e., — 4^ + v(x), and its adjoint. We establish 
: H ■ the equivalence of this description with the mathematicians' definition of spectral singularities 

$_i ■ for the potential v(x) = z_5(x + a) + z + 5(x — a), where z± and a are respectively complex 

& '. 

and real parameters and 5(x) is the Dirac delta-function. We offer a through analysis of the 
spectral properties of this potential and determine the regions in the space of the coupling 
constants z± where it admits bound states and spectral singularities. In particular, we find 
an explicit bound on the size of certain regions in which the Hamiltonian is quasi-Hermitian 
and examine the consequences of imposing 'PT-symmetry. 
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1 Introduction 



The use of non-Hermitian Hamiltonians in theoretical physics has a long history. It extends from 
early attempts to construct divergence- free relativistic quantum field theories [TJ, to more practical 
and successful applications in nuclear and atomic physics [2] and particularly quantum optics 
[31 H]. During the past ten years there has been a renewed interest in the study of a special class of 
non-Hermitian Hamiltonians that possess a real spectrum. The best-known examples are the VT- 
symmetric Hamiltonians [5] such as p 2 + ix 3 . These belong to the wider class of pseudo-Hermitian 
Hamiltonians H whose adjoint H* is given by 

H^ = r]H7]-\ (1) 

for some Hermitian invertible operator t], 0, [7J. What makes pseudo-Hermitian Hamiltonians 
interesting is that they are Hermitian with respect to a possibly indefinite inner product^, namely 

<V>, := <■!»?•>, (2) 

where (-|-) is the inner product of the Hilbert space in which all the relevant operators act. Most 
of the recent work on the subject is concentrated on a particular class of pseudo-Hermitian Hamil- 
tonians, called quasi- Hermitian [8], that satisfy (JTJ for some positive-definite (metric) operator rj. 
In this case, (J2J) is a positive-definite inner product, and H becomes Hermitian provided that we 
define the physical Hilbert space of the system using the inner product (•, -) v , [9j[T0]. This allows for 
formulating the pseudo-Hermitian representation of quantum mechanics in which PT-symmetric 
as well as non-'PT-symmetric quasi- Hermitian Hamiltonians can be employed to describe unitary 
quantum systems [10]. The techniques developed in this framework have so far found interesting 
applications in relativistic quantum mechanics [11], quantum cosmology [12], quantum field theory 
[13] . bound state scattering [H], and classical electrodynamics [T3]. But these developments do not 
undermine the importance of the requirement that the observables of a unitary quantum system 
must be Hermitian with respect to the inner product of the physical Hilbert space, [TO] . 

Among the properties of Hermitian operators that make them indispensable in quantum me- 
chanics is their diagonalizability. For Hermitian and more generally normal operators (those com- 
muting with their adjoint), diagonalizability is equivalent to the existence of an orthonormal basis 
consisting of the eigenvectors of the operator. This is more commonly referred to as completeness. 
For a non-normal (and hence non-Hermitian) operator H, diagonalizability of H means the exis- 
tence of a basis B^ consisting of (scattering and bound state) eigenf unctions of the adjoint operator 
iff that is biorthonormal to some basis B consisting of the eigenfunctions of H, i.e., B and B^ form 
1 This means that {-,H-) V = (H-, ■),,. 
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a biorthonormal system of the Hilbert space, [10]. For brevity we shall call such a biorthonormal 
system a biorthonormal eigensystem for H. 

Diagonalizability is a weaker condition than Hermiticity, but diagonalizable operators with a 
real and discrete spectrum can be related to Hermitian operators via similarity transformations [7J . 
This in turn implies that they are quasi-Hermitian [5]. The situation is more complicated when the 
spectrum is continuous. A serious difficulty is the emergence of spectral singularities that conflict 
with the diagonalizability of the operator in question [16]. The aim of this paper is to elucidate 
the mechanism by which spectral singularities obstruct the construction of a biorthonormal eigen- 
system for the operator. We shall achieve this aim by obtaining a quantitative measure of lack 
of a biorthonormal eigensystem and comparing the latter with the mathematical condition for the 
presence of spectral singularities that is based on the behavior of the Jost functions. In order to 
clarify the meaning and consequences of spectral singularities we shall offer a detailed investigation 
of the spectral properties of the Hamiltonian operators of the form 



where (± are complex coupling constants, a is a real parameter, and 5(x) stands for the Dirac 
delta-function. 

An alternative mechanism that can make a non-Hermitian operator non-diagonalizable is the 
emergence of exceptional points. These correspond to degeneracies where both eigenvalues and 
eigenvectors coalesce [17]. Exceptional points have various physical implications [31 IH [18]. But 
they must not be confused with spectral singularities. Unlike exceptional points that can be present 
for operators with a discrete spectrum (in particular matrix Hamiltonians) , spectral singularities 
are exclusive features of certain operators having a continuous part in their spectrum. As we 
will see in Sections 2 and 3, for an operator having a spectral singularity we can still define two 
linearly independent (scattering) eigenf unctions for each eigenvalue, nevertheless it is impossible 
to construct a biorthonormal eigensystem for the operator. To the best of our knowledge, physical 
meaning of spectral singularities and their possible applications have not been previously studied. 
This is the subject of Ref. [19] where the results of the present paper have been used to develop a 
physical interpretation for spectral singularities. 

Ref. [20] uses the mathematical theory of spectral singularities developed in [2U E2] to empha- 
size their relevance to the recent attempts at using complex scattering potentials to define unitary 
quantum systems. The results of [20] are, however, confined to potentials defined on the half-line 
x > 0, where the Hamiltonian operator acts in the Hilbert space of square-integrable functions 
ip : [0, oo) — > C satisfying the boundary condition ip(0) = 0. Furthermore, due to the nature of the 
concrete potentials studied in [20] , it has not been possible to construct bases of the corresponding 
Hamiltonian and its adjoint and show by explicit calculation how the presence of a spectral sin- 




(3) 
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gularity obstructs the existence of a biorthonormal eigensystem. This is quite essential, because 
for the cases that the spectrum is real, the availability of a biorthonormal eigensystem is a neces- 
sary condition for the existence of an associated metric operator and the quasi-Hermiticity of the 
Hamiltonian, [TO] . 

The only thoroughly studied example of a complex scattering potential that is defined on the 
whole real line and can lead to spectral singularities is the single-delta-function potential with a 
complex coupling [23]. The Hamiltonian operator is given by fl3]) with a = £_ = 0. It develops 
a spectral singularity if and only if the coupling constant (£+) is imaginary. In particular, for the 
cases that the real part of ( + is positive, the bound states are also lacking and the Hamiltonian is 
quasi-Hermitian. The complex single-delta-function potentials provide a class of manifestly non- 
PT-symmetric Hamiltonians with a continuous spectrum that happen to be quasi-Hermitian. An 
advantage of considering complex double-delta-function potentials is that their space of coupling 
constants has a subspace, given by ( + = d, where the Hamiltonian is PT- invariant. Therefore, 
these potentials provide an opportunity to investigate the significance of PT-symmetry [23] . 

The spectral properties of the PT-symmetric double- and multiple-delta-function potentials 
have been studied in [25l [26l [271 [28] . The results are, however, confined to the determination of 
the (scattering and bound-state) spectrum of these potentials, and no attempt has been made to 
decide if these potentials lead to spectral singularities. 

In the present paper we will try to obtain a map of the space C 2 = M 4 of the coupling constants 
C± that specifies the regions corresponding to the existence of bound states and spectral singular- 
ities. We will in particular investigate the intersection of these regions with the two-dimensional 
PT-symmetric subspace: £_ = C+- The following is an outline of the results we report in this 
paper. In Section 2, we obtain an explicit quantitative measure of the existence of biorthonormal 
eigensystems and compare the latter with the known condition of the presence of spectral singular- 
ities. Here we also provide a useful characterization of spectral singularities and bound states for 
complex scattering potentials. Section 3 treats the spectral properties of the double-delta-function 
potentials. It consists of four subsections in which we obtain the regions in the space of coupling 
constants where spectral singularities and bound states exist, find their location in the spectrum of 
the operator, and determine a lower bound on the size of certain regions in C 2 where the operator 
(j3~l) is quasi-Hermitian. Section 4 presents our concluding remarks. 

2 Spectral Singularities 

Consider a complex-valued potential v : R — > C depending on a set of complex coupling constants 
Zi,z 2 , - ■ ■ ,Zd such that v* is obtained by complex-conjugating the coupling constants in the expres- 
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sion for v. Suppose that v decays rapidhcl as \x\ — > oo and that the spectrum of the corresponding 
Hamiltonian operator, 
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d 

H=-— + v(x), xeR, (4) 

dor 

is the set of nonnegative real numbersjj Let ipl k (x) denote the (generalized) eigenfunctions of H, 
i.e., linearly-independent bounded solutions of 

Htfjx) = k 2 r ak (x), (5) 



where k G M + and a G {1,2} are respectively the spectral and degeneracyj labels and z : = 

By definition, H is diagonalizable, if ipl k {x) together with a set of (generalized) eigenfunctions 
4> z ak (x) of H* form a complete biorthonormal system {ip* k , 4>l k }, i.e., they satisfy 

2 



2 /"OO 

K> = ^ *(* - g), X / ^ l^iX^Sbl = 1, (6) 

a=l ■ /o 



where (-|-) is the usual L 2 -inner product. The biorthonormality relations fl6]) imply the spectral 
representation of H, 

H = J2 dkkV ak )(cf>t k \, (7) 



0=1 

as well as the eigenfunction expansion: 



2 ^oo 

f( x ) = ^ / rf /c /a* $*0e), 
0=1 - 70 



where / : R — > C is a test function and 

/«* == (Al/>- (9) 

Because iff = — -J^ + t>(x)*, are the eigenfunctions of . This in turn means that cf) z ak 
must be a linear combination of he., there are J a b G C satisfying 



6=1 



(10) 



2 For the purpose of the present paper, we assume that as \x\ — > oo we have |v(x)| < exp[— eyjaij] for some e e R + . 
As far as the general properties related with spectral singularities, all the results hold for the less rapidly decaying 

potentials that satisfy dx (1 + < oo. See [29l I30j . 

3 Wc shall consider the more general case that the spectrum involves eigenvalues (with square-integrable eigen- 
functions) at the end of this section. 

4 The spectrum is necessarily doubly-degenerate. 
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In view of §6§, there must exist K ab £ C such that 

(€k\r bq ) = K ab 5(k-q). (11) 

Furthermore, if we respectively denote by /, J and K the 2x2 identity matrix and the matrices 
with entries J ab and K ab , we find J*K = I. In particular, K must be an invertible matrix and 
J a b = Kab*- We can write this relation in the form 

K* 

Jab = d^(fp (12) 
where K is the transpose of the matrix of cofactors of K. It satisfies 

= K5(k-q). 

k \n q ) J 

We can use (1121 and (ITUl) to express ([9]) as 



/«* = 3i^^^«*<^£l/>- ( 13 ) 



According to this equation if det(iT) = 0, the eigenfunction expansion (jBJ) breaks down; the eigen- 
functions do not form a complete set; and H is not diagonalizable. We identify this situation 
with the presence of spectral singularities: 



Spectral singularities are points k 2 of the continuous spectrum of H where det(K) = 0. (14) 



Because of (complex) analyticity property of the eigenfunctions ipl k1 det(i^) is an analytic 
function of k. Therefore the (real) zeros of det(K) are isolated points forming a countable subset 
of the real line. Moreover, because v is a bounded function decaying rapidly away from zero, 
the eigenfunctions tend to plane waves as k becomes large. This shows that det(K) does not 
have arbitrarily large zeros (for fixed z). As a result, the zeros of det(-fT) are not only isolated 
but actually finite in number. In other words, depending on the values of the coupling constants 
zi,Z2, - ■ ■ ,Zd, det(K) may have no (real) zeros in which case spectral singularities do not arise and 
H is diagonalizable, or a finite number of (non- vanishing real) zeros K\, K2, ■ ■ ■ , in which case 



^ are spectral singularities and H is not diagonalizable. 



In general the space of coupling constants can be divided into two regions, namely the singular 
region where H has spectral singularities and the regular region where it is diagonalizable. 
In mathematics literature a spectral singularity is defined as follows: 

Definition 1: An element E* of the (continuous) spectrum of H is called a spectral singularity 
if the integral kernel of the resolvent operator: (H — E)" 1 , i.e., the Green's function (x\(H — 
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E) 1 \y) , is an unbounded function in every small open neighborhood of E+, but E* is not an 
eigenvalue of H with a square-integrable eigenfunction 



There is a rather general theory of spectral singularities for the differential operators of the form 
P| where the spectral singularities are characterized as the real zeros of certain analytic functions 
[HI EH [321 [331 E2l [29]. For the case that the operator acts in L 2 (R), this is the Wronskian, 

W[^,^ k+ ] := - i/>' k _(x)if> k+ (x) = ^-(0)^ + (0) - Vfc-(0)^fc+(0), (15) 

of the Jost solutions ipk± of the eigenvalue equation Hip = k 2 ip. These are defined in terms of their 
asymptotic behavior 

ipk-(x) -> e~ lhx for x ->■ -oo, i>k+(x) ->■ e tkx for x ->■ oo. (16) 

More specifically, we have 



Spectral singularities are the real (non-vanishing) zeros of W[ipk-, fpk+\- (^) 



This description of spectral singularities seems to differ from the one given in (JHj). In Section 3, we 
demonstrate the equivalence of the two descriptions for the double-delta-function potential by ex- 
plicit calculations. The following calculation shows how the description (I14p relates to Definition 1. 
Using p, ©, flU and flT2J, we have 

{X \(H - E) -'\y) = ± r* j&m^L = ± r dk mmm: 

o=l ^° a,b=l ^° 

K ab ^ k (x)^f k ( y y 



" poo 

E / 

- [.-1 -/o 



det(iT)(A; 2 — E) 



o,b=l 

In the remainder of this section we provide a useful characterization of the spectral singularities 
and bound states. 

Because \v{x)\ decays rapidly as \x\ —> oo, solutions of (jSJ) have the asymptotic behavior: 

^ ka {x) -> A ± e ife * + £ ±e - ifcE for a; -)• ±oo, (18) 

where A± and 5± are possibly fc-dependent complex coefficients. If we denote the coefficients A± 
and B± for the Jost solutions ipk± by A± and B±, we can express (fT5|) as 

= fil = 1, = fi+ = 0. (19) 

Next, we let M = (M a t,) be the possibly fc-dependent 2x2 transfer matrix [M] satisfying 



5 Therefore spectral singularities are certain poles of (x\(H — E) 1 \y). 
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and use this relation and Eqs. (ITS]) and ffT9|) to obtain 



A + = 5+ = -^^7, AT = M 12 , B7 = M 22 . (21) 

~ detM' detM' + 12 ' + 11 y 1 



Inserting these equations in ( fTSl) . we find 

^ fc _(x) M 12 e ifcx + M 22 e _ifcx for x ->■ oo, (22) 
^ fc+ (x) -»■ -— ^ for x -oo, (23) 

Because according to Abel's theorem [35], the Wronskian of solutions of (jSJ) is independent of x, 
we can use the asymptotic formulas for the Jost solutions to compute their Wronskian. We use 
Eqs. (|16p . 022 p and (]23p . to perform this calculation first for x — > oo and then for x — > — oo. This 
gives 

W^fc-^fc+l = 2ikM 22 (k), (24) 
W ^~^ = "detW' (25) 



where we have made the /c-dependence of M 22 and M explicit. A direct consequence of (l24j) and 

d25D is 

detM(fc) = 1. (26) 
More importantly, we have the following characterization of spectral singularities that follows from 

mm and m 



Spectral singularities are given by k 2 where k is a (non-vanishing) real zero of M 22 (k). (27) 



Finally, consider the more general case that the Hamiltonian operator fl4]) has, in addition to a 
continuous spectrum corresponding to k G M + , a possibly complex discrete spectrum. The latter 
corresponds to the square-integrable solutions of (JHJ) that represent bound states. It is not difficult 
to show that the spectral label corresponding to these bound states are also zeros of M 22 (k), but 
unlike the zeros associated with the spectral singularities these must have a positive imaginary 
part. In other words, we have the following characterization of the bound states. 



Bound state energies are given by k 2 where k is a zero of M 22 (k) with $s(k) > 0. (2 



3 The Double-Delta-Function Potential 



3.1 Eigenfunctions for Scattering and Bound States 



Consider the time-independent Schrodinger equation 

h 2 d 2 



2m dx 2 



+ C+5(x-a) + C-5(x + a) 



V> = Erp. (29) 



Let £ be an arbitrary length scale and introduce the dimensionless quantities 

2m£(± x a 2m£ 2 E . . 

Then ( 1291) takes the form 

- ^" + [z+6{x - a) + z_S(x + a)]il> = Eip. (31) 
We can write solutions of fl3"TT) as 

(ip~(x) for x < —a, 
ij°(x) for \x\<a, (32) 
^ + (x) for x > a, 
r(x) = A u e lkx + B u e~ lkx , ue{-,0,+}, (33) 

where k := V^E and without of loss of generality we require that the principal argument of k belongs 

tO [0,7T). 

To determine the matching conditions at x = dba, we demand that ip be continuous, i.e., 

ip-(-a) = ifj°(-a), ifj°(a) = ip + (a). (34) 

Furthermore, we integrate both sides of (13T|) over the intervals [=Fa — e, =Fa + e] and take the limit 
e — >■ in the resulting formulas to find 

^"'(-°) - ^°'(- a ) + z-ip°(-a) = 0, tp '(a) -tp + '(a) + z + tp°(a) = 0. (35) 

Introducing 

w ± := — , (36) 

and inserting fl32|) and fl33l) in (pl|) and fl35|) yield the desired matching conditions that we can 
write in the form 

\ / V _ w e -2iafc j_ _ w _ J y Bq J ' y B + J ~ \ w + e 2iak l + w + ) \ B J ' 

(37) 

In light of these relations, the matrix M satisfying ( 120]) reads 



M 



1 — W- — w + + (1 — e 4iak )w-W + 2iw-iu + sin(2a/c) — w-e 2tak — w + e 2tak 
-2iw-W + sin(2a/c) + W-.e~ 2tak + w + e 2mk 1 + w_ + w + + (1 — e 4iak )w-W + 

It is easy to check that indeed det(M) = 1. 
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Next, we let z stand for (z-, z + ) and use ip± k and ip% k ^° denote the eigenfunctions obtained by 
setting A = ^tt)" 1 / 2 , B = and A = 0,B = (2tt)^ 1 / 2 respectively. Then, 

^f fc (x) = (2tt)~ 1/2 x { e ikx for Id < a, (39) 



i + 


^ikx _ w e ~ik(x+2a) 


for 


x < —a 






for 


\x\ < a. 


l-w+] 


j e ikx _|_ w ^ e -ik(x-2a) 


for 


x > a, 


1 — wJ 


j e ~ikx _|_ w e ik(x+2a) 


for 


x < —a 






for 


\x\ < a. 


i+w + : 




for 


x > a. 



i>l k {x) = (2tt)- 1/2 x { e ~ ikx for \x\ < a, (40) 



We can construct a set of eigenfunctions of H* by taking z± to z± or w± to — w± in these relations. 
They are given by 

i/jf k (x) = (2tt)- 1/2 x { e ikx for |x| < a, (41) 



1 - w*_ 


)e lkx + w*_e 


—ik(x+2a) 


for 


x < —a 








for 


\x\ < a, 


l+w* + 


)e lkx — w* + e 


—ik(x—2a) 


for 


x > a, 


l+w*_ 


)e~ lkx — w*_ 


gifc(x+2a) 


for 


x < —a 








for 


\x\ < a, 


1-wl 


\ e -ikx _|_ w * 


gifc(x— 2a) 


for 


x > a. 



ip^ k (x) = (2tt)- 1/2 x { e ~ ikx for |x| < a, (42) 



3.2 Characterization of Spectral Singularities 

In this subsection we use (fT4|) to determine the spectral singularities of the double-delta-function 
potential. This requires computes (^ffcl^bg) f° r an a > ^ e {1, 2}. Using (I3"9"j) and (HOI and the 
identities 



le 



e^ x dx = tt5(h) + , / e^ x dx = n6(fi) 

we find 

where K = (Kij) is a 2 x 2 matrix with entries 



S(k-q)K, (43) 



z_ + Zj 



K u = K 22 = l-wt-wl = l+ ~ 4fc2 + , (44) 

K 12 = w_(l - w_)e 2iak - w + (l + w + )e~ 2iak 

= (4k 2 )- 1 [iz-{2k - %zJ)e 2iak - iz + {2k + iz + )e~ 2mk ] , (45) 

K 21 = -w^(l + w„)e~ 2iak + w+(l - w+)e 2mk 

= (4k 2 )- 1 [-iz_{2k + iz_)e~ 2iak + iz + {2k - iz + )e 2iak ] . (46) 
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In the "PT-symmetric case, where z + 

An = 
K 12 = 
A 21 = 



K 



22 



--: z 7^ 0, K is a real matrix, and 
1 ■ 



2k 2 



(2k 



[*(2fc + 



jz e 



-2iak 



(2A; 2 )" 1 3 [z(-2fc + zz)e 2iafc ] 



(47) 
(48) 
(49) 



The fact that if is not generally diagonal shows that {V'afc? ^Pfq} * s n °t a biorthonormal system. 
To construct the basis biorthonormal to {ipl k } we transform ipf k according to 



6=1 



and fix the coefficients J a b by demanding that {ip% k , <fi z ak \ be a biorthonormal system. As we ex- 
plained in Section 2, in terms of K this condition takes the form 5 ab = J^=i ^Acb- Therefore, a 
basis biorthonormal to {"0^1 ex ists provided that the matrix K is invertible, and the matrix J of 
coefficients J a t, has the form J = K~ u . 

The nonzero real values of k for which if is a singular matrix give the spectral singularities of 
H . These are the non- vanishing real zeros of det(if) that we can obtain using 



det(if) = 1 + 



z 2 _ + z 2 
4k 2 



+ 



z 2 z 2 , 



8k 4 



+ 



z_z + 
2k 2 



Ak 2 



cos(4a/c) 



2_ + Z^ 

2k~ 



sin(4a/c) 



0. (50) 



If either Z- 



and z := z + or z + = and z := Z-, this equation reduces to 

;2 



z 
4k 2 



0. 



Therefore, for pure imaginary z there is a spectral singularity located at k — ±iz/2 = \z\/2. This 
agrees with the results for the single delta-function potential 
For the 'PT-symmetric case (z + = z*_ =: z), we have 



1 ~ 4F ' cos ^ 4a/c ) 



2fc 2 8fc 4 2k 2 
In particular if z is purely imaginary, i.e., z = ia for some a G 



sin(4afc) 



0. 



(51) 



det(AT) = cos 2 (2a k) 



1 



2k 2 ) 



sin 2 (2aA;). 



(52) 



Therefore, det(if) = iff cos(2a/c) = and k = \a\/y/2. This implies that there is a spectral 
singularity for k = \o\j\[2 = \z\/y2 iff a takes one of the following values 



cr, 



7r(2n + 1] 
2V2 a 



neZ. 



(53) 
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In summary, for the case that z + = —z_ =: z is purely imaginary, H has a single spectral singularity 
given by 

rr 2 [(On 4- 1 Wl 2 

(54) 



= o$ = [(2^+ 1)tt 
2 



4a 

if z = ia n for some n G Z. Otherwise it does not have any spectral singularities. 

Next, consider the general "PT-symmetric case where z + = z*_ =: z, and z need not be purely 
imaginary. In this more general case, we rewrite (IBTj) in the form 



where 



It is easy to compute 



f(z,a,k) :-- 

K[f(z,a,k)] = 
$*[f(z,a,h)] = 



det(K) = \f(z,a,k)f 



k| 2 sin(2afc) +e _ 2iak 



2k 2 



$t{z) 



\z\ 
2k 2 



l\ sin(2a£;) + cos(2a£;), 



cos(2a/c) H — sin(2a/c) 



(55) 
(56) 

(57) 
(58) 



In view of (|55l) . det(if) = iff a, &)] = $s[f(z,a,k)] = 0. Imposing £y[/(z, a, fc)] = 0, we 

have 

(59) 



cos(2a/;;) = — sin(2a/c), 



k 

which in particular implies sin(2afc) ^ 0. Moreover, cos(2a/;;) = iff dt(z) = 0. In light of 
sin(2a£;) ^ and ([59}, 9£[/(z, a, fc)] = gives 



- 9?(z) 2 + 3(z) 2 = 2fc 2 . 

This implies that if > ^(z)!, there is no spectral singularity. 

Next, we solve (1601 for k to obtain 



(60) 



k 



■$t(z) 2 + %(z) 2 



(61) 



and express (1391 as 

&( z ) = -fccot(2afc). (62) 

Inserting flBTT) in fl62|) yields a necessary and sufficient condition for the existence of a spectral 
singularity, namely 

2»(z) tan (aV2[-3?(z) 2 + 3(z) 2 ]) + ^[-^(z) 2 + 3(z) 2 ] = 0. (63) 
Introducing the variables 

r:=2a9f?(z), s := 2a3f(z), t := a^/2[-^(z) 2 + %z) 2 }, (64) 



12 



-40 - 20 20 40 

Figure 1: Curves in the r-s plane giving the location of the spectral singularities 
for the general PT-symmetric double-delta-function potential. The dashed 
lines are the asymptotes s = ±r. The intersection of the curves with the 
s-axis corresponds to the spectral singularities given by Eq. (1541) . 



we can express fl63|) in the form 

r = —t cot t, (65) 

and establish 



s = ±tVl + csc 2 t. (66) 

Figure [H shows a plot of the parametric curve defined by ( 1651) and ( |66l) . It consists of an infinite 
set of disjoint open curves with asymptotes s = ±r in the r-s plane. The points on these curves 
correspond to the values of the coupling constant z for which a spectral singularity appears. 

Next, consider the general not necessarily PT-symmetric case. Generalizing our treatment of 
the "PT-symmetric case, we use (15"Uj) to factorize det(i^) as 

det(K) = f_(z-,z + ,a,k)f+(z-,z+,a,k), (67) 

where 

f ± {z-,z + , a, k) := ^ sin(2a£;) + e ±2 * ak ± i) , u := z_z + , v := Z —^ ± - (68) 

Therefore det(K) = if and only if at least one of f±(z~, z+, a, k) vanishes. 
Let us abbreviate f-(z_, z + ,a,k) as f(k), i.e., set 

f(k) := ^ sm(2ak) + e~^ k g - i) . (69) 
13 



Then it is easy to see that f+(z_, z + , a, k) = —f(—k). Therefore, the positive zeros of f+(z-, z + , a, k) 
are identical with the absolute- value of the negative zeros of f_(z_,z + ,a,k). In other words, the 
spectral singularities are given by positive and negative real zeros of f(k). Another interesting 
property of f(k) is that it satisfies 

f(k) = -ie~ 2iak [I + + w+ + w^w + {l - e 4aik )} = -ie~ 2iak M 22 {k) , (70) 

where M 2 2{k) is the the entry of the matrix M of (138]) with the row and column labels 2. According 
to (1701) the spectral singularities are the non- vanishing real zeros of M 2 2{k). This establishes the 
equivalence of (TT4]) and (1T7I) for the double-delta-function potentials. 

In order to characterize the real zeros of f(k), we set the real and imaginary parts of the 
right-hand side of ( 169]) equal to zero. This gives 

1 + H r-^ sm(2ak) + — 



j sin(2a£;) + ^ ) cos(2ak) = 0, 
^ ) \in(2afc)+f^-l)cos(2afc) = 0. (72) 



2k 2 k J y \ k 

Because sin(2afc) and cos(2a/c) cannot vanish simultaneously, these equations hold provided that 
the matrix of coefficients of sin(2a/c) and cos(2a/c) is singular. Equating the determinant of this 
matrix to zero and simplifying the resulting equation, we find 

g(k) := k 3 - 2%(v)k 2 + + |^ k + \ [SR( W )^( W ) - SR(^)S(w)] = 0. (73) 

Because g is a real cubic polynomial, it always has at least one real root k. If k ^ 0, E* = k 2 
is a spectral singularity. Expressing k as a function of u and v and inserting it in say ( 1721) we 
find a sufficient condition on the coupling constants z± for the existence of a spectral singularity. 
Repeating this for all the roots of g (for the cases that (1731) has other nonzero real solutions) we 
obtain a complete characterization of the spectral singularities. They lie on a three-dimensional 
surface S embedded in the four-dimensional space (C 2 ) of the coupling constants (z-, z + ). Figured] 
is a graphical demonstration of the intersection of S with the plane z + = that represents the 
"PT-symmetric region of C 2 . In the following we examine some non-"PT-symmetric regions of C 2 
and their intersection with S. 

1. Consider the plane Hi in C 2 defined by z + = — z*_ =: z where u = —\z\ 2 and v = i^s(z). Then 
Eqs. UZIJ and flZ2J take the form 



M(z)\ 2 , (Q( Z ) ^ 



k J V k 



sm(2ak) = = ( — ^ - 1 ) cos(2aA;). 



These are satisfied if and only if sin(2a/c) = and k = \^s(z)\. Therefore, we have a spectral 
singularity located at E* = k 2 = ^s(z) 2 if and only if 

^s(z) = nir, for some n G Z — {0}. (74) 
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This shows that Hi intersects S along equidistant lines parallel to the 9fJ(,z)-axis in III. 

2. Consider the case that both z + and z_ are purely imaginary. This also defines a plane in C 2 
that we denote by IT2. In this case, we can express z± as 

iy± 

Z± =: , 

a 

where y± are nonzero real numbers. In terms of y±, Eqs. 07ip and (]72|) take the form 

y+ + y- «\ /r> 1 \ n fy+ + y- y+y- S\ . M ,\ 

1 cos(2afc) =0=1 : — — 1 sm(2afc). 



2ak J V 2 « 2fc2 

There are two ways to satisfy these equations. Either 

y+ + y 



- 1 = sin(2a£;) = 0, (75) 
2ak 

or 

2a 2 k 2 - ak (y + + y_) + y + y_ = cos(2ak) = 0. (76) 
We consider these two cases separately. 

If (1751) holds, E* = k 2 with k := (y + + y^)/2a is a spectral singularity provided that 

y+ + y- = ™> ™eZ-{0}. (77) 
This defines a set of equidistance parallel lines in n 2 along which we have spectral singularities. 
If (|76p holds, k = n n where K n := (2n + l)ir/ (4a) for all n G Z, and 



2/+ = a«; n 



?/_ — 2<2K n \ (2n + 1)7T 



2y_ - (2n + 1)tt 



(7* 



y_ — a«; n / 2 L%- — (2^ + I) 71 " 

This equation gives the location of another set of spectral singularities, namely E* = k 2 , in 
the plane II2. 

Figure [2] shows the curves in n 2 along which a spectral singularity arises, i.e., II 2 H 5. 

3. Consider the plane II3 in C 2 corresponding to z + = —z_ =: z. Then v = and w = — z 2 . 
In particular, 5s(u) = — 2^R,(z)5s(z) . We can confine our attention to the subcase: $s(u) 7^ 0, 
because for $s(u) = either ^t(z) = 0, in which case z± are purely imaginary and the results 
of case 2 apply, or 5s(z) = 0, in which case the potential is real and there are no spectral 
singularities. 

In view of v = and (1731) . g(k) = k(k 2 — 9ft(w)/2). Therefore, k does not have a real zero and 
there is no spectral singularities, if 9ft(w) < 0. For 9ft(w) > 0, there is a spectral singularity at 
E„ = 4 = &(u)/2, where 

K± - (T 9 ) 
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Figure 2: Curves in the y+-y~ plane (Il 2 ) along which one has a spectral 
singularity for purely imaginary couplings. There are spectral singularities 
along the y_- and y + -a,xes. The dashed line (y + = —y~) represents the VT- 
symmetric double-delta function potential with purely imaginary couplings. 
The intersection of these lines with the full curves correspond to the spectral 
singularities given by Eq. fl54|) . 



Inserting this equation in (172|) gives 

%(u) = ± Sft(u) cot (ay/m(uty . (80) 

Introducing the parameter t := a-\/29ft(u), we can use flHUl) to obtain the following para- 
metric equations for the r := a$t(z) and s := aQ(z) values that correspond to the spectral 
singularities. 

\r(t)\ = lV\csct\ - 1, \s(t)\ 



2 a/ | sin 1 1 — sin 2 1 

Figure |3] shows the graph of the parametric curves defined by ( 18TT) . They form the intersection 
of the plane II3 with the singular region S of C 2 . 

4. Consider the case that z± — (1 + is±)/a with s± e ffi. arbitrary. This corresponds to another 
plane in C 2 that we denote by 1I4. Introducing 

S- + S+ l + s_s+ 

s ■= 2 ' * := 2 ' ^ ' 

we have 

1 + is 1-t + is 

v = , u = . (83) 

a a z 
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Figure 3: Curves in the r-s plane along which spectral singularities occur for 
the coupling constants with opposite sign. The origin (s = r = 0) does not 
actually lie on these curves. The intersection of the curves with the s-axis 
corresponds to the spectral singularities given by Eq. (1541) . The dashed lines 
are the lines s = ±r. 



Inserting these in (173]) yields 

a 3 g(k) = (ak - s)(a 2 k 2 - ask + t) = 0. (84) 
Therefore, we need to consider the following two possibilities. 



4. a) k = s/a = (s_ + s + )/(2a). In this case ( !72|) is satisfied automatically while (17T]) 
yields 

t = scot(2s) + 1. (85) 
We can use (182|) and (|85l) to express s± in terms of s. This gives 

s_ = st v / s 2 + l-2(scot(2s) + l), (86) 
s+ = 2s- S - = s± ^s 2 + 1 - 2(scot(2s) + 1). (87) 

4.b) k 7^ s/a. Then according to (1841) . 

t = ask- a 2 k 2 . (88) 
Furthermore, both (jTTl) and (j72l) become 

tan(2aA;) + aA; = (89) 
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Figure 4: Curves in the s_-s + plane along which the spectral singularities occur 
for the coupling constants of the form z± = l+is±. The solid (red) and dashed 
(blue) curves correspond to the spectral singularities with k = (s + + s_)/2 
(Case 4. a) and k ^ (s + + s_)/2 (Case 4.b), respectively. Also shown (in green) 
is the unit disc: s 2 _ + s 2 + < 1, where there are no spectral singularities. 

This equation has a countably infinite set of real solutions K n that can be easily obtained 
numerically. Substituting K n for k in (1581 and using ( 152|) . we find 



Figure H] shows the parametric plot of the curves in the plane corresponding to the 

spectral singularities for both cases 4. a and 4.b with a — 1. As seen from this figure, there 
are no spectral singularity in the unit disc defined by s 2 _ + s 2 + < 1. 

3.3 Location of the Spectral Singularities and the Bound States 

As we noted in Section 2, the spectral singularities are given by the real zeros of M22(k) while the 
bound states correspond to zeroes of M 2 2{k) with positive imaginary part. For the double-delta- 
function potential, we can write M 22 (A;) = in the following more compact form. 



s =F \Js 2 + 1- 2(asn n - a 2 ^), 



(90) 
(91) 



2s — s_ = s ± a/s 2 + 1 — 2{asK n — a 2 K 2 ). 



(£-3_)(Jt-3 + ) 



(92) 



where we have used ( 1381) and introduced 



2ma(± 



& := 2%ak. 



3± := az± = 
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In particular, the spectral singularities are give by 

K := (93) 
where .ft is a nonzero solution of fl92|) lying on the imaginary axis in the complex .ft-plane, i.e., 

Re £ := {w G C | 9ft(w) = ^ 10 }, 

whereas the bound state "energies" are given by ( I93p for solutions .ft of ( 19 2 j) lying to the left of this 
axis, i.e., 

^en_:={weC 9ft(w) < o }. 

For both spectral singularities and bound states, we have 9ft ( .ft) < which implies \e m \ < 1. 
Taking the modulus of both sides of (!92|) . we find 

|£-3_||ii-3 + | < |3+||3_|. (94) 

This is violated for any .ft fulfilling 

>\l-\ and |^-3+| >|3+|- (95) 
Therefore, the solutions of (192]) with 9ft(.ft) < must belong to the union of the discs 

£>± := {Re C | |i?-3±| < |3±| }• 

This provides an upper bound on the size of the region in the complex .ft-plane where bound state 
energies and spectral singularities are located, namely 

%:=(n„uf)n(fl + uD_). 

Here we have used the index 3 := 3+) to emphasize the 3±-dependence of 1Z^. Figure |5] illustrates 
the discs D± and the region 7Z^ for a generic choice of 3± and also for the case that %± are real and 
positive. It is easy to see that in the latter case is empty and there are no spectral singularities 
or bound states. 

Let D a and D a be the disc and half-disc defined by 

D a := {Re C I \R\ < a }, D CT := {.ft G C | |.ft| < a, 9ft(j?) < 0. }, (96) 

where a is the largest of 2|3±|, i.e., 

a:=2max(|3_|,|3 + |). (97) 

Then, .ft G So- is a weaker necessary condition for the existence of bound states and spectral 
singularities. This is simply because D± C D a . See Figure EI a)- 
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Figure 5: (a) Discs D± and -D CT for generic values of %±. The gray area with the 
origin excluded is the region 1Z? where the bound states and spectral singularities 
are located (if any), (b) D± for 3 ± £ M + . In this case is empty. 



Because the spectral singularities and bound states are given by the zeros of 

:= (A - 3 _)(* - H) - 3-3+ e 2 \ (98) 

which is an entire (everywhere complex-analytic) function, and these zeros are contained in D a 
which is a compact subset of the complex ^-plane, we can determine the location of the spectral 
singularities and bound states using the following well-known result of complex analysis: 

Theorem 1: Let C be a counterclockwise oriented contour bounding a compact and simply- 
connected region R in complex plane and h : C — > C be a function that is analytic on an 
open subset containing R. Then h has a finite number of zeros in R. Moreover, if none of 
these zeros lie on C, the contour integral 

^=^#4 (99) 
2m J c h(w) 

gives the sum of orders of zeros of h contained in R. In particular, if all of these zeros are 
simple (of order 1), ric gives their number, [361 §10]. 

A proper use of this theorem requires a careful analysis of the order of zeros of F^. It is not 
difficult to show that the zeros of can at most be of order three. Moreover, A is a third order 
zero of Ff if and only if & = and 

— 1 ± i 1 — 1 =F i , 

3-=— 3+ = ^ = ^- (1»°> 
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This does not correspond to a spectral singularity or a bound state. Ff has a second order zero £2 
if and only if 



2 3 _ 3+ e l+3-+3 + ±Vl+(3-3 + ) 2 = 1 ± v /l + (3_-3+) 2 , (101) 

1 



£2 = 2 [1 + 3- + 3+ ± V / l + (3--3+) 2 J • (102) 

Requiring that 9ft (£ 2 ) < 0, we can use (110 II) to show that |1 ± a/1 + (Jl — 3+) 2 l < 2|3_3 + |. There- 
fore, there is no spectral singularity or bound state associated with a second order zero of Ff, if for 
both choices of the sign, 



1 ± v/l + (3_- 3+ ) 2 > 2|3_ 3+ |. (103) 

This inequality in turn implies the following sufficient condition for the lack of spectral singularities 
and bound states associated with a second order zero of Ff. 

13-3+1(13-3+1 -1)< l3 '~ 3+ ' 2 - (104) 

In particular, such bound states or spectral singularities are forbidden if |3_3 + | < 1. 

Next, we return to the idea of using Theorem 1 for locating the spectral singularities and bound 
states. For this purpose we can use the contours C(p, 9) and c±(p±) depicted in Figure[6]to compute 

1 / fJW 



n(p-,p + ) := n-(p-) + n+(p+), (106) 

1 r Fi(&) 

N(p,9) := — j> -jfr^dSt, (107) 

where p, p± G [e, a], e <C 1, and 9 G [e,n — e]. In the generic case where Ff has no second order 
zeros, n±(p) and N(p, 6) give the number of zeros of Ff enclosed by c±(p) and C(p, 9), respectively. 
Therefore, plotting n(p_, p+) and N(p, 9) as functions of p± and (p, 9), we can locate all the spectral 
singularities and bound states of the double-delta-function potential for given coupling constants 
z±. In particular, for e — > 0, n tot := n(a,a) and N tot := N{a,n — e) respectively give the total 
number of spectral singularities and bound states, except for the cases that for some imaginary 
^ both ^ and — ^ are zeros of Ft. In the latter case, R and — ^ give rise to the same spectral 
singularity, and one must account for the corresponding double counting in n tot . 

Note that locating spectral singularities is most conveniently carried out using (1731) . In the 
absence of an analogous equation giving the k values for the bound states, we use N(p, 9) to locate 
the latter. Figure [7J shows the density plots of N(p, 9) for the PT-symmetric case 3+ = — 8 ± 3i 
and the non-PT-symmetric case 3_ = —8 + 3i, 3+ = —4 — 2i. These resemble the phase diagrams 
of statistical mechanics where the critical points correspond to the bound states. As we expect, for 
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A. 




c+(p+) 



c-(P-) 



(a) 



(b) 



Figure 6: (a) c±(p±) are rectangular contours of width e < 1 and hight p± < a; 
(b) C(p,9) is the boundary of the region lying between circular arcs of side 
length e C 1 and p < a and opening angle 6 G [e, re — e] . 



the PT-symmetric case the location of these points is symmetric about the real axis in the complex 
^-plane. 

Figure [8] shows the graphs of N(p, it — e) for the PT-symmetric case 3± = — 1 ± 8i and the non- 
PT-symmetric case 3_ = — 2 + 7i, 3+ = — 4 — 5i. These show the distance between the bound states 
from the origin. For the PT-symmetric case the bound states are created in complex-conjugate 
pairs with the same distance from the origin. This explains the fact that the number of bound 
states changes in increments of 2. This is clearly not the case for the non-PT-symmetric case. For 
both of the above choices of the coupling constants, a < 17. Therefore, the maximum value of each 
curve gives the total number of bound states for the corresponding system. 

Figure shows a contour plot of N tot := N(a, ir — e) for 3 ± = 1 + is± and e = 10~ 6 as functions 
of s± G R. Although the real part of the coupling constants are positive and equal, for large enough 
values of their imaginary part the system develops bound states. This is in contrast to the single 
delta-function potential where there are no bound states for coupling constants with a positive real 
part. We also see that in the PT-symmetric case s + = — s_, which corresponds to the depicted 
diagonal line, the number of bound states change in increments of 2. This is consistent with the 
fact that these are produced in complex-conjugate pairs. 

3.4 Real Bound States and Quasi-Hermiticity 

An important feature of the graphical demonstration of the location of spectral singularities and 
bound states in the complex ^-plane is that for the cases that 9ft(3±) > and |Q : (3±)| are sufficiently 
small, the system does not have any spectral singularities or bound states. Figures H] and [9] provide 
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Figure 7: Density Plots of N(p, 9) for the PT-symmetric case 3± = —8 ± 3i (on the 
left) and the non-"PT-symmetric case 3_ = —8 + 3i, 3+ = —4 — 2i (on the right) in 
the complex .ft-plane. & r and 84 mark the real and imaginary axes. As the color 
changes from the lightest to the darkest N(p, 9) takes values 0,1, and 2, respectively. 
The critical points marked by black spots are the ^-values corresponding to bound 
states. They are symmetric about the ^.-axis for the PT-symmetric case. 

a clear demonstration of this phenomenon for the case that 9?(3±) = 1. 

The presence of spectral singularities is an obstruction to the quasi-Hermiticity of the Hamil- 
tonian operator. This is also true for the bound states unless they happen to have real energies 
(eigenvalues). We will refer to these bound states as "real bound states." It is not difficult to see 
that generic bound states are not real. In this subsection, we shall first derive analytic expressions 
for the existence and location of real bound states and then for fixed and positive values of 9ft(3±) 
we establish the existence of a positive lower bound on the size of a region in the Q : (3„)-S(3 + ) plane 
where the system is free of both the spectral singularities and bound states. This is a region where 
the Hamiltonian operator is quasi-Hermitian. It is in this region that we can employ the machinery 
of pseudo-Hermitian quantum mechanics [HI [10] to construct an associated positive-definite metric 
operator and use the Hamiltonian operator to define a unitary quantum system. 

Excluding the case of a single-delta-function potential [21] where 3+3. = 0, we can express (|92|) 
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Figure 8: Plots of N(p,ir — .01) (top figure) and N(cr,9) (bottom figure) for the 
■PT-symmetric system defined by i± = — 1 ± 8i (the solid curves) and the non-PT- 
symmetric system defined by 3_ = — 2 + 7i and 3+ = — 4 — 5i (the dashed curves). 
For the PT-symmetric model iV(p, ir — .01) changes in increments of 2 while N(a, 9) 
is symmetric with respect to the 9 = ir/2 line. 



as 

(HS- 1 )-* (108) 

According to fl93l) . the real bound states are given by real and negative solutions of (11081) . For these 
solutions the right-hand side of f ll08p is real, positive, and less than 1. Equating the left-hand side 
with its complex-conjugate yields 

3(3-3+)* = l3-| 2 3(3 + ) + |3+| 2 ^)- (109) 

In order to explore the consequences of this equation we introduce the notation 

r± := 3*(3±), s ± := 3( 3± ), (110) 

and consider the following cases separately. 

(i) 3(3_3 + ) = 0: In this case, 

r_s + + r + s_=0, ^— + ^±- = 0. (Ill) 

l3-r l3+r 

Therefore, either both s± vanish and the potential is real or both s± are nonzero. In the 
latter case, ( lllll) implies 3_ =3+- This is the 'PT-symmetric case for which ( J92l reduces to 

|^-3+| = |3+|^. (112) 
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Figure 9: Contour plot of the number iV tot of bound states located in the region: 
7r/2 + e < arg(^) < 37r/2 — e for $± = l + s±i and e = 10 -6 . As the color changes from 
the lightest to the darkest N tot take values 0,1,2,3,4, respectively. The diagonal 
line s + = — s_ corresponds to the PT-symmetric region along which the number of 
bound states change in increments of 2. 



Because e R < 1, this equation cannot be satisfied, if 9ft(3+) > 0. This is consistent with 
the results of [28]. Furthermore, for the non-PT cases with real 3_3 + , such as 3_ = —3+ or 
imaginary %± with 3_ 7^ 3+, there is no real bound states. 

(ii) S(3_3 + ) 7^ 0: In this case, we can write (j!09|> as 

|3-| 2 ^(3 + ) + l3 + | 2 ^(3-) ril ,x 

Substituting this equation in (]9"2]) gives a rather complicated relation between 3_ and 3+. 
This relation together with the requirement that the right-hand side of (j!13p be negative 
provide the necessary and sufficient condition for the existence of real bound states for non- 
"PT-symmetric cases. We have implemented this condition to address the existence of real 
bound states for the special cases where 3+ = }-e lu := 3 e tu ^ 2 with v E [0,2n). Figure [TUI 
shows the curves in the complex 3-plane along which real bound states exist for various values 
of v. It is important to note that all these curves are finite in length. Therefore, there are no 
real bound states for sufficiently large values of |3|. 

Next, we wish to show the existence of regions in the space of the coupling constants %± where 
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Figure 10: Curves in the complex 3-plane along which real bound sates 
exist for j ± = %e ±iu/2 , v = 7i(2n - l)/20, and n G {1, 2, • • • , 10}. The 
numbers attached to each curve segment is the corresponding value of 
n. r and s respectively mark the 9? (3)- and 3(3)-axes. Note that all the 
curves have finite length. 

there is no spectral singularities or bound states. Our main tools are the following basic theorems 
of real and complex analysis. 

Theorem 2: Let n G Z + , D be a compact subset of K n with its standard topology, and 
Lp : M. n — > M. be a function that is continuous on D. Then {<p(x) \ x G D } has both a 
minimum and a maximum, [23 §8]. 

Theorem 3 (Maximum Modulus Theorem): Let C be a contour bounding a compact 
and simply-connected subset R of the complex plane and h : C — > C be a function that is 
analytic on an open subset containing R. Then {|/i(w)| | w G R } attains its maximum on 
C, [381 §111-1]. 

First, we use Theorem 3 to prove the following preliminary results. 
Lemma 1: Let D p denote the following half disc of radius p G M + : 

D p := { R G C I |£| < p, < } , 
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and L : C — > C be the function defined by 

f 1 - e m 

L{&):=\ M f ° r (114) 

[ -2 for £ = 0. 

Then \L\ attains its maximum value on D p at .ft = 0, i.e., 2 = |L(0)| is the maximum of 
A p := {\L(&)\ | R G D p } for all p G M+. 

Proof: First, consider the case p < 1. Then 2) p C Di, which implies A p C ^/t^. Therefore, 
the maximum of A p is less than or equal to that of A\. This shows that it is sufficient to 
prove the lemma for the case p > 1. Because L is an entire function and D p is compact, 
according to Theorem 3, A p has a maximum that is located on the boundary of D p . This 
is the union of the closed line segment £ p := {iy \ y G [— p, p] } and the open semicircle 
C p := {ip e tip | (f G (0, 7r) }. The maximum of A p is the largest of the values taken by \L\ on 
£ p and C p . We will show that these values are bounded from above by 2. Because G D p 
and |-^(0)| = 2, this is sufficient to prove the lemma. In the following we consider the values 
of \L\ on t p and C p separately. 

— For all & G £ p , we can write & = iy for some y G [— p, p}. Inserting & = iy in (11141) and 
computing the modulus of both sides of the resulting expression yields 

. . . . 2 sin y 

\L(&)\ = — — - < 2. (115) 

— For all & G C p , we can write &= ip e lif for some ip G (0, tt). Because p > 1 and sin tp > 0, 
f[TT4D implies 

= |l-exp(2zp e^)| ^ 1 + exp(2zp = j + e -2 P sin V < 2 _ (n6) 
P 

This together with (11151) proves the lemma for p > 1. As we explained above this 
establishes the statement of the lemma also for the case p < 1. □ 

Lemma 2: Suppose that r± > 0. Then ^ = is a first order zero of the function defined 
by flHED- 

Proof: Recall that the zeros of are at most of order three and ^(0) = 0. Therefore, 
it is sufficient to show that 8. = is not a second or third order zero of F^. Assume (by 
contradiction) that R = is a second order zero of Fp Then Fl(0) = 0, i.e., 3_+3_+2j_3 + = 0. 
Equivalent ly, we have 

r + + r_ + 2(r_r + — s_s + ) = 0, s + + s_ + 2(r_s + + r + s_) = 0. 
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Solving the second of these for s_ and inserting the result in the first, we find 

2(1 + 2r„) S 2 

r + + r_ + 2r_r + H = 0. 

1 + 2r + 

But this equation cannot be satisfied for r± > 0. This shows that the above assumption is 
false and ^ = is not a second order zero of Fp Next, we recall that ^ = is a third order 
zero of Fj* if and only if fllOOp hold. But these conflict with the condition r± > 0. Hence 
.ft = is not a third order zero of Fp □ 

Next, we use Theorem 2 and Lemmas 1 and 2 to prove the following desired result. 

Theorem 4: Suppose that r± > and \s±\ < r max := max(r_,r + ). Then there is a positive 
upper bound Bp on \s±\ such that for all s± satisfying \s±\ < Bp, the Hamiltonian ([3]) does 
not have any spectral singularities or bound states^ 



Proof: Recall that spectral singularities and bound states are zeros .ft of F$ with 9? (.ft) < 
and that they belong to D a , where a := 2 max(|3_|, |3+|). The latter is a subset of the half-disc 



&|<V8r max , »(£)<()}, 



because in view of r± < r max and \s±\ < r max , we have a < \/8r m3X . According to Lemma 2, 
.ft = is a first order zero of This implies that the function Gf : C — > C defined by 

, . f iT^Oft) for & 0, . . 

3 \ Fi(0) for £=0, 

is an entire function and Gf(O) ^ 0. Furthermore, the spectral singularities and bound states 
of the Hamiltonian ([3]) correspond to the zeros &q of lying in D. Another important 
observation is that Gy has no zeros .ft with 9? (.ft) < 0, because if existed these zeros would 
have corresponded to the spectral singularities or bound states of the Hamiltonian ([3]) with 
real and positive coupling constants ($± G But as we argued above this Hamiltonian 

does not have any spectral singularities or bound states. This observation establishes the fact 
that 

GV(£)^0, for all .ft G £. (118) 

Because GV is an entire function, |GV| is continuous on D which is a compact subset of C = M 2 . 
In view of Theorem 2, this implies that the set {|GV(.ft)| | .ft G 3) } has a minimum m?, i.e., 
there is -ft min G T) such m? = \Gp(R min )\. Because .fto,.ft m in G T) and (11181) holds, we have 

< |GV(£ mi n)| = m? < \G?{8v)\. (119) 



3 Here and in what follows r := (r_,r_|_). 
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Next, we introduce : C — > C as the function defined by 

J ? (£):=G ? (£)-G r -(£). (120) 

Because Gj(Ro) = 0, we have 

|J ? (^o)| = |G r ,(i?o)|. (121) 
Furthermore, in view of (fTTTj) . (jI20j) . (jlHjl . and the fact that .ft ^ 0, 

J r (^o) = + s+) + [ - + z(r_ S+ + r + s_)]L(^o). (122) 

This implies 



|^o)| < + \s+\ + [\s-\\s + \ + |r_||s+| + |r+||s_|J|L(^o)| 

< 2 (3r max + 1) s max , (123) 

where s max := max(|s_|, \s+\) and we have used the triangular inequality, the condition 
|s±| < r max; and |L(^o)| < 2 that follows from Lemma 1. 

If we combine fl 123[) with (11211) and (11191) . we obtain 



771 ~* 

< 2<s~r+T) £ — ■ < 124 > 

This inequality is violated for the values of s± for which 

< - < 125 > 

Therefore, for the cases that \s±\ < B? the existence of .fto leads to a contradiction; such a .fto 
cannot exist; and there are no spectral singularities or bound states. □ 

The upper bound B? given in (11251) involves the minimum m? of |Gy| on the half-Disc 2). 
Because Gy is a nowhere- zero analytic function on 2), l/Gy is also analytic on D. Hence, according 
to Theorem 3, 1/|GV| attains its maximum on the boundary of D. It is not difficult to see 
that rrtf = 1/Mp. Therefore, in practice, for given values of r±, we can obtain by exploring the 
values of |Gy| on the boundary of D. 

We can identify the boundary of D with the graph IV of the parameterized curve: 

7r -(t) := 2zr max [{2t + l)6(-t) + e i7Tt Q{t)] , £ G [-1, 1], (126) 

where 9 is the unit step function: 6(0) := 1/2 and 6(t) := (1 + t/\t\)/2 for t ^ 0. Figure [TT] shows 
the graphs of \G?{pi?{t))\ and |L(7^(t))| for the case r± = 1 that is considered in Figures 0] and [HI 
In this case, r max = 1 and D is the half-disc of radius 2 lying in Il_. As seen from the graph of 
values of \L\, it attains its maximum at t — —0.5 (corresponding to .ft = 0) and its maximum value 
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Figure 11: Plots of |C^(7^(t)) | (the full curve) and \L(jf(t))\ (the dashed curve) 
as a function of t G [—1,1]. £_ ~ —0.949 and t + ps —.051 are the minimum 
points of |G^(7^(t))| corresponding to R± ~ ±1 .795i. These give the minimum value 
m r w 1.906. The maximum value of |L(7^(t))| is 2 that is attained at to = —0.5 
corresponding to A = 0. 

is 2. This is consistent with the statement of Lemma 1. The minimum points of |GV| are located 
at t — —0.949 and t = —.051. These correspond to ^ m i n ~ ±1.795z where \G?\ takes its minimum 
value: w 1.906. According to (11251) . this gives = ~ 0.238. Therefore, for $± = I ± is± 
with \s±\ < 0.238 there should be no spectral singularities or bound states. This is in complete 
agreement with the graphical data depicted in Figures H] and [U the disc with center s± = and 
radius 0.238 lies in the region with no spectral singularities or bound states. 

4 Concluding Remarks 

In this article we provided an explicit demonstration of how spectral singularities obstruct the 
existence of a biorthonormal eigensystem and render the Hamiltonian non-diagonalizable. We 
achieved this by obtaining a characterization of spectral singularities in terms of the M 22 entry of 
the matrix M of Eq. (1201) . In particular we showed that while bound states are zeros of M 2 2{k) 
with Q(k) > 0, the spectral singularities are the real zeros of M22(/c). It is not difficult to infer 
from this observation that, similarly to the bound states, the spectral singularities are linked with 
singularities of the scattering matrix [19J. However, unlike the bound states, they lie on the real 
axis in the complex /c-plane. This in turn suggests interpreting spectral singularities as resonances 
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having a vanishing width. Ref. [19] provides a thorough description of this interpretation and its 
physical implications. 

We established the utility of our general results by providing a thorough analysis of the spectral 
properties of a two-parameter family of complex point interactions. We obtained various results on 
the nature and location of the bound states and spectral singularities for this family and proved 
the existence of regions in the space of coupling constants where both bound states and spectral 
singularities are lacking and the Hamiltonian is quasi-Hermitian. 

Throughout our study we examined the consequences of imposing "PT-symmetry which cor- 
responds to restricting the coupling constants to a complex plane in the space C 2 of coupling 
constants. This revealed a previously unnoticed fact that PT-symmetric double-delta-function 
potential can involve spectral singularities. 

The results of this paper may be extended to complex point interactions corresponding to three 
or larger number of delta function potentials. Another line of research is to try to compute a metric 
operator r] + and the corresponding equivalent Hermitian Hamiltonian h and the pseudo-Hermitian 
position and momentum operators X and P for the double-delta-function potential whenever the 
Hamiltonian is quasi-Hermitian. Theorem 4 provides the mathematical basis for a perturbative 
calculation of r] + , h, X, and P. We plan to report the results of this calculation in a forthcoming 
publication. 
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